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I. INTRODUCTION 

The Haystack antenna pointing system makes use of a coordinate 
conversion program (COCON). This program accepts celestial inertial 
coordinates (right ascension, declination, distance, and time) referenced to 
the equator and equinox of date of a particular point (e.g., star, planet, or 
artificial earth satellite) and produces the pointing angles and range to that 
point from the Haystack installation. Range rate is also calculated using the 
celestial coordinates and their rates. 

A second program called RADEC is used in the system for the purpose 
of deriving right ascension and declination from the range and pointing angles. 
This program drives the right ascension and declination display lights on the 


master control console. 


ia MA THEMA TICAL RELA TIONSHIPS 
A. Basic Relationships 


An ellipsoidal earth is assumed with an equatorial radius of A = 1. 


(See Fig. 1). The flattening is taken as f = 357° The polar radius is 
B =A(1—f). Given the geodetic latitude, O.> of a site, the geocentric 
x 
latitude O is found. 
es meee cea ear. (1) 
gE Eg 


* The following notation is used in this report. A subscript E will denote 
reference to the earth anda subscript s will denote the object under 
observation. 


The radius of the earth to the site (at 0 height) is 


A sec On 


Re = ae (2) 
E eens tan” O, | 


The perpendicular distance to the radar plane, located at a height h, above 


E 


the surface of the earth, from the geocenter is 


te = ty + Ry cos O, — 6, ) (3) 


The distance ap to the vertical from the geocenter is 


d,. = Rye sin G — 6, ) (4) 
B. Inertial Celestial Coordinates to Radar Cartesian 
Coordinates 


Given the celestial inertial coordinates, the geocentric cartesian 
coordinates of the object are found. The Z axis goes through the North 
Pole, the X axis through the vernal equinox, andthe Y axis makes a 


right handed orthogonal system. 


X = R cos 46 cos @ 
S S S S 
= R cos 4 sin@ 
S S S S 
Z = R sin 6 (5) 
S S S 


where R. is the distance from the geocenter to the object, 6. is the 
declination, and a is the right ascension from equinox. 

The origin of the radar cartesian coordinate system is the radar site. 
The Z axis is vertical, the Y axis points north and the X axis points east. 


The radar cartesian coordinates are found by the matrix relationship 


etl, 
x 


x. i 0 0 —sin or cos Sa 2 0 
Y.. = {0 sin 9 cos o. - |—cos Sh —sin Qe. O; - Y + ds. (6) 
| “ai O —cos 9. sin 9. 0 0 l \*s Th 


where sb is the local sidereal hour angle of the radar site and is found hy 


om = D(ST, —ST,) + ST. + Ap (7) 


where D is the fractional part of the observation time expressed in 


days, ST, and ST. are the apparent sidereal hour angles of 


1 
Greenwich for the midnights (00 hrs. GMT) succeeding and preced- 
ing the time of observation, and Ae is the longitude of the site 


(positive east). 


The expansion of Eq. (6) leads to the following expressions 


X 


_ Ro cos oe sin G = 
. R,| sin 6. cos 9. —sin 9. cos 6. cos (a, —&,) | + oe 


Z R,| cos oe cos 9. cos (a, 2, ) + sin 6. sin 1 | ae (8) 


Y 


Ce Radar Range From Celestial Coordinates 


Taking the square root of the sum of the squares, the range R is 


found from Eq. (8): 


ReeeR [1+ 
S 


R 
S 


(9) 


Although Eq. (9) appears somewhat cumbersome, the values in the paren- 


thesis involving Do da and ry. are stored as constants. Note that in the 


case of distant stars R, oo and hence R-> R.- 


D. Azimuth From Celestial Coordinates 


In order to find azimuth, first an angle AZI is found using Eq. (8) 
x 
AZI = tan™} 


| 


Ne 


r 


-] cos 6. sin (a, —2, ) 


ti 
ctr 
+) 


(10) 
qe 
sin on cos 0. — cos 6. sin 0. cos (a, 2, ) + a 
Ss 


The azimuth (A), which is measured from north, is then determined using 
the following rules: 


ie X and Mes are both positive, A lies in the first quadrant 


andA =AdZI 


i xX. is positive and _- is negative, A lies in the fourth 


quadrant and A = 180° —AZI 


it x. and a are both negative, A lies in the third quadrant 
and A = 180° + AZI 
if x is negative and Y. is positive, A lies in the second 


quadrant and A = 360° ~AZI 


If es Ss the numerator of (10) goes to zero giving rise to the 


following special cases: 
O 
= = 
when 6. Dn A=0 


6. < Do A = 180 


Oe D> object directly overhead. 


An examination of the denominator of (10) leads to the following 


Special cases: 


as 6. = 90°. the denominator — » and the arc tan of AZI -0. 


The conclusion is that 
if 6. = 90°, A=0°, 


if 6. = —90°, A= eel 


If the denominator = 0, then for 


if 
iN @) 
© 
0 


a, 7 & , A 
a < 2 , A= 270, 


a, =p, the point is directly overheadandA is arbitrarily set to 0°. 


E. Elevation From Celestial Coordinates 


The elevation (E) is found using the results of Eqs. (8) and (9) in 


1 “R 


k= sin a (11) 


As the numerator goes to zero the elevation angle goes to zero. If 


the numerator is negative the object is below the radar horizon, 


F, Range Rate From Celestial Coordinates 


The expression for computing range rate is found simply by 


differentiating Eq. (9). 


R K, sin 6 K, cos 6, cos (a4, —%), . 
; l S PA S 


+K, sin 6_ cos (a, —2,)} b. + K, cos a sin (o. —Q) (x. —2,,) | (r2) 


where Ky and K, 


of the celestial coordinates are provided as input to the coordinate conversion 


are constants (defined in Section III Bla). The derivatives 


program while the derivative of the local sidereal hour angle of the site is 


derived from the apparent sidereal hour angles. 


G. Right Ascension and Declination From Radar Coordinates 


Solving (6) for the geocentric cartesian coordinates 


X —sin aa —sin B.. cos Sogn cos 0. cos a x. 
< = cos aa —sin § sin, cos o. sin Y-dp (13) 
Z| 0 cos 0. sin Wo Zi "E| 


In this case Ks Yo and 2. are found from the radar pointing angles 


x = Rcos E sinA 
= Rcos E cos A 
Z = R sinE (14) 


where R is radar range, E is elevation, and A is azimuth. 
Expanding Eq. (13) and using the relationships in (14) the celestial 


cartesian coordinates are found: 


di. sin an + Tp Cos Ny 


| | | = 
AL=R cos 2, {sin E cos 9.—cosE cos A sin, + SS 
—cos E sinA sinQ | = RX! 
E Ss 
d_ sin9.. + r_cos @ 
. ; : E E E E 
Y. =R Sin $2, | sin E cos 9. — cos E cos A sin Be + +=} 
+ cos E sinA cos Q2 | = RY' 
E Ss 
d.. cos 9. — rw sin 9 
Z -=Ricos Ecos A cos 9. + sinE sin §, — E E E EE] 
S E E R | 
= ' ] 
= RZ° (15) 


Right Ascension is found by 


eee ee 
O = tan S 
ve 
Ss 


[xy + Gry” 
S a 
In the event R~o, Eqs. (16) and (17) still hold since 


limit R 


Roc : 
II. THE COORDINATE CONVERSION PROGRAM (COCON) 
A, COCON Inputs 


There are four types of inputs to COCON 


1. Common Storage Parameters 


(16) 


(17) 


Site parameters and physical constants are found in common 


storage. The following quantities are used: 
Site geodetic latitude in degrees with B20 
Site longitude in degrees with B20 


Site height in feet with BO 


Equatorial Earth Radius in nautical miles with B17 


¥ Conventionally a B of 20, or B20, means the binary point is to the right of 


bit 20. 


2. Greenwich Sidereal Time 


The first file of the ephemeris tape* contains the apparent side- 
real time of Greenwich. This file is comprised of two records. The first 
file gives the day number of the first entry of the second file, the epoch 
date being January 1.0, 1964. The second record contains one word for 
each day corresponding to 00 hours GMT. This word gives the apparent 
sidereal time of Greenwich in radians with a B of 26. 

a Console Writer 

If the sidereal times for the required days are not found on the 
ephemeris tape, or if a tape error is encountered, the experimenter is 
given the option of having the program try to read the tape again or of having 
it compute** the sidereal time. 

4, Inertial Celestial Coordinates and Rates 


The arguments are: 


SRA Right ascension in revolutions B27 
SDEC Declination in revolutions B27 
RADIUS Distance - if in earth radii the distance is 


stored with B22 

- if in astronomical units the dis - 
tance is complemented and 
stored with B24 


- if mo, a zero is stored 


#See Group Report 1964-41, "Haystack Pointing System: Ephemeris Tape 
Program", D.M. Hafford, 25 September 1964 


##For the year 1964, the sidereal time of Greenwich in radians for the ue 


day is given by ST. = i(. 01720279) + 1. 7193827. 


RADOT Time rate of change of right ascension in 


radians /sec B37 


DECDOT Time rate of change of declination in 


radians/sec B37 


RA DIOSDOT Time rate of change of radius in nautical 
miles/sec B24 


CONVERTIME GMT of observation in days B28 
Be COCON Program Details 
ie Initialization 


The initialization section performs the following functions: 
a. Computes Tr and do from Eqs. (3) and (4) 


Computes: K, = 2(cos 9.4, — sin irre) 
K, =e sin $d, — cos Inte ) 
K, = (rg +4, ) 
2, = (ST,,,—ST, )/86400 


KIDOP = Ky in nautical miles for range rate 


computation 
K2DOP = K, innm 


b. Read from the ephemeris tape or compute 
ST, STi ST. where i is the day of the experiment 
(January 1.0 = 1) 
c. Set local sidereal hour angle using CONVERTIME. 
ae Working Section 
The working section of COCON is a straight forward computation 
program using Eqs. (9), (10), (11) and (12) to find range, azimuth, elevation 


and range rate. 
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C. COCON Outputs 


COCON outputs left in common storage are: 


GEOCENLAT geocentric latitude in degrees B20 
SIDER TIME local sidereal hour angle in radians B26 
YRTRAN de. from Eq. (4) 

ZRTRAN Tp from Eq. (3) 

AZIM azimuth in revolutions B27 

ELEV elevation in revolutions B27 

RANGE in radar units with BO 

TRUERANGE += ER radii B22 


complement inA.U. B24 


0 = STAR (arbitrarily infinite) 


IV. THE RADEC PROGRAM 
ee RA DEC Inputs 
We Common Storage 


RADEC makes use of the site parameters and some constants 
computed by COCON. These are: 
YRTRAN, ZRTRAN, aot Am? FRAMESIZE 
a Typewriter 
The experimenter has the option of selecting a set of coordinates 
from which right ascension and declination are derived. There are two dis- 


tinct classes of coordinates available as input to RADEC: 


11 


a. Celestial as found in common storage 

RA and DEC in revolutions B27 

RA and DEC both with scan in revolutions B27 

ley Azimuth and elevation 

. Actual the actual azimuth and elevation 
of the antenna - an infinite range 
is used, 

Command the azimuth and elevation com- 
mands being sent to the antenna - 
infinite range. 

Uncorrected AZ- azimuth, elevation, and true 

jap 8 range as output from COCON. 

. AZ-ELtSCAN azimuth and elevation as modi- 
fied by the scan program, and 
true range. 

Corrected AZ - azimuth and elevation with scan 

EL as corrected for site characteris - 


tics, and true range. 


B. RADEC Program Details 


ie Initialization 
The initialization section of this program computes two constants, 
sets the time delay, and sets the switch for selecting the required input 


coordinates. 


iz 


The computed constants are Ky and K, and are used in Eq. (15) 


~ 


= dp. sin Oe + Tp COs Oo 


~ 
M" 


3 op cos Oe —rp, sin Or 


Since pointing angles are computed ahead of real time to allow 
for interpolation, the local sidereal hour angle as stored by COCON does 
not correspond to the time for the actual radar coordinates. In order to 
compensate for this fact, the amount of earth's rotation in 2 framesizes 
is subtracted from SIDERTIME. 

The selection of input coordinates to RADEC is at the discretion 
of the experimenter. When the pointing system is cranked up, RADEC 
selects actual radar coordinates. The experimenter may, via the typewriter, 
modify RADEC to select any one of the above mentioned 7 inputs. 

2. Working Section 

The working section of RADEC operates once every second or 
once aframe. If the system is cycling normally (a frame every two seconds), 
RADEC's working section operates every second. The first entry is from 
the Master Control Program, while the second entry is effected via the one 
second interrupt. If the system is cycling in a high speed planning mode 
(no computer time buffering between frames), the interrupt portion of 


RADEC is disabled. 
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The rest of the working section is a straight forward computation 
routine. If the inputs are radar coordinates, the expressions found in IIG 
are used to generate right ascension and declination. If the inputs are in 


celestial coordinates, coordinate conversion is not required. 


Ss RADEC Outputs 

There are two RADEC outputs. The first output uses channel 5* to 
drive the display lights. This is accomplished by an external-function com- 
mand. Figure 2 shows the format of the words going out on channel 5. The 
right ascension is displayed in time (hours, minutes, seconds) whereas the 
declination is in degrees (degrees, minutes, seconds). The local hour angle 
is sub channel 6 and has the same format as the right ascension. (The local 
hour angle is, SIDERTIME - right ascension. ) 

The second output stores right ascension and declination in ASTRORA 
and ASTRODEC for use by the radiometer program. These registers are 


the images of the words going out on channel 5. 


*J.E. Gillis, A.F. Dockrey and S.B. Russell, to be published. 
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Fig. 1. Ellipsoidal Earth 
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